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Abstract 

In this paper, we calculate the transmission and reflection amplitudes of wave func- 
tions for different potentials such as the delta function, the rectangular barrier, the 
Eckart potential, and the Hulthen potential. We describe the relationship between 
these amplitudes and compute the reflection resonances between each potential. We 
describe the transmission and reflection probabilities using the WKB formula and 
compare the results with ones obtained from matching the boundary conditions. Fur- 



thermore, we use a two by two transfer matrix to calculate a rigorous bound on the 
transmission and reflection probabilities. 
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1 Introduction 

One-dimensional quantum problems widely appear in any textbook on quantum mechanics 
[D[3l3jlS[g[6j[3|8j|g[Tnj[^ Although they are simple in mathematics 

and clear in physics, they lead to important and new results [T7J [18j [19j I2Q1 [2TJ [22j 123] . 
In physics, there are a number of applications of one-dimensional quantum problems. For 
example, in acoustics, one might be interested in the propagation of sound waves down a 
long pipe, while in electromagnetism, one might be interested in the physics of wave-guides 
[23]. 

In this paper, we study one-dimensional quantum problems for four selected potentials, 
the delta function, the rectangular barrier, the Eckart potential, and the Hulthen potential 
|24j . We will calculate exact and approximate transmission and reflection probabilities by 
matching boundary conditions and by using the WKB method respectively, and compare 
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them to each other. We will also use a 2 x 2 transfer matrix to obtain a rigorous bound 
on these probabilities. Moreover, we will derive resonances of transmission and reflection 
probabilities. 



2 Conventions 



We are interested in solving the time-independent Schrodinger equation [U 121 El lU El El d 

is ia am nu aai us], 

~^ + V{x)]^{x) = E^x). (1) 
2m dx z 

We concentrate on a potential which is asymptotically constant, that is V(x) — > as 
x — > ±oo. For E > V(x) for all x, the asymptotic wave functions are given by 

exp (-ik +oc x) 



t- 



i/j(x) — > 



+oo 



exp (—ik_ 00 x) exp (ik^^x) 
h r- 



where wave numbers are defined by 



k 2 



Vk- 



2m (E - V ±OG ) 

h 2 



if x — > +oo 
if x — > — oo 



(2) 



(3) 



We identify t and r as transmission and reflection amplitudes respectively. Therefore, the 
conservation of flux leads to the condition 

T + R = l, (4) 

where T = |t| 2 and R = \r\ 2 are transmission and reflection probabilities respectively. A 
position at which a transmission or reflection probability is in unity, it is thus referred to as 
the occurrence of resonance. 

Moreover, we can obtain a rigorous bound on the transmission probabilities by using a 
2x2 transfer matrix. It gives 



-^transmission ^ Secll ( 



X-2 



XI 



k ~ 



k 2 {x) 



dx 



where 



k 2 {x) 



2m [E - V(x)} 

h 2 



(5) 



(6) 



and k(x) — > k if x is outside the interval (xi,x 2 ). 

On the other hand, for E < V(x) in the range x\ < x < Xi we can obtain the transmission 
probability for any potential by using a useful technique called the WKB method. The result 
is given by 



exp 



(7) 
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3 One-dimensional problems 



In this section, we review main results in some one- dimensional problems. Four types of 
potentials are selected to discuss in this paper. 

3.1 A delta function potential 

A delta function potential takes the form 

V(x) = a5(x), (8) 
where a is a positive constant and the delta function is defined by 

We are interested in a scattering state E > 0. We define 

, 2mE , , ma , , 

k =-fp- and ko = -tf- (10) 

The transmission and reflection amplitudes for this potential are given by (see 011]) 

t = k and r = tk ° ■ (11) 
k — ikq k — iko 

Therefore, the transmission and reflection probabilities are 

T = WT¥ and R= ¥T¥- (12) 

Note that 

T + R = l. (13) 

The transmission and reflection probabilities varying with k are shown in Fig. [U We 
have seen that the transmission probability tends to unify as k goes to infinity. We say, 
however, that this potential has no transmission resonances. On the other hand, the reflection 
resonances occur at k = 0. Moreover, from the plotting if the potential is strongest, the total 
reflection happens. That is when the potential is stronger, penetration of a particle or wave 
through the potential is more difficult. 

3.2 Rectangular barrier potential 

The rectangular barrier potential has the form (see [21 U\) 

| 0, if otherwise ' 

The shape of the rectangular barrier is shown in Fig. [2j We are interested in two cases 
E > V > and V > E > 0. 
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Figure 1: Plotting of transmission and reflection probabilities varying as k for the delta 
function potential with (a) ko = 1, (b) k = 2, (c) k = 10, (d) k = 100 and (e) k = 1000. 



4 



3.2.1 Case I: E > V > 

We define 

2 2mE 2 2m(E-V ) 2 2mV o ,22 / 1K x 

fc = g = ^ , and A; = = k - q . (15) 

The transmission and reflection amplitudes for this potential are given by (see [TBI 126] ) 

4fcgexp(2z/ca) . . 

(16) 



(k + g) 2 exp(2iga) — (k — g) 2 exp(— 2zga) 

and 



2i(Ar — g 2 ) sin(2ga) exp(2i£;a) 



(A; + g) 2 exp(2iga) — (k — q) 2 exp(— 2zga) 
Therefore, the transmission and reflection probabilities are 

r 4fc V and ^ fc 4 sin 2 (2ga) 

4A; 2 g 2 + ^sin 2 (2ga) 4Pg 2 + k% sin 2 (2ga) ' 1 J 

Note that 

T + i? = l. (19) 

The transmission and reflection probabilities varying with q are shown in Fig. [3] and 
HJ Fig. |3] compares the effects of barrier heights Vq of the potential on the probabilities 
when the width a of the potential is fixed. It has been found that the higher the barrier of 
the potential is, the more the number of reflection resonances. This is similar to the case 
of the delta function potential: penetration of a particle or wave through the potential is 
hard to occur when the barrier height of the potential is large. Fig. H] compares the effects 
of the widths of the potential on the probabilities when the barrier height of the potential 
is fixed. The results are that the reflection resonance can occur when the width increases. 
Analytically, the transmission resonances occur at 

9 =2? (20) 

where n — 1,2, 3, while the reflection resonance is at k = for this potential. 

Moreover, we can obtain a lower bound of the transmission probability by using a 2 x 2 
transfer matrix 

r£sech2 (v*?)' (21) 

This rigorous bound is compared with the exact solution in Fig. [5] It is found that the lower 
bound approaches the exact solution when a particle is at high energy. 
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3.2.2 Case II: V > E > 

We define 



, 2m (V - E) 



Q = — ^ ( 22 ) 

The transmission and reflection amplitudes for this potential are given by 

2iQke~ 2ika , s 

(k 2 - Q 2 ) sinh(2Qa) + 2ikQ cosh(2Qa) 1 ; 



and 

(k 2 + Q 2 ) smh(2Qa)e- 2ika 



(k 2 - Q 2 ) sinh(2Qa) + 2ikQ cosh(2Qa) 
Therefore, the transmission and reflection probabilities are 



(24) 



Ak 2 Q 2 

{k 2 - Q 2 ) z sinh 2 (2Qa) + 4k 2 Q 2 cosh 2 (2Qa) 



T = 3 - ^ - (25) 



R = ^ W 1 ) (26) 



and 

kl sinh 2 (2Qa) 
{k 2 - Q 2 f sinh 2 (2Qa) +4PQ 2 cosh 2 (2Qa) 

We have seen that in this case T ^ 0. It follows that a particle can penetrate the barrier 
from one side to the other although the potential energy of the particle exceeds its total 
energy, which does not appear in classical physics. This is called tunneling. Note that in 
this case 

T + R^l. (27) 

By the WKB method, we obtain 

T w = exp(-AQa). (28) 

In Fig. [6j the result from the WKB method is compared with the exact result. It is found 
that the larger the height of the potential, the more accurate the WKB approximation 
becomes. 

3.3 An Eckart potential 

The Eckart potential takes the form 

v{x) = —2— + —2— tanh U + (29) 

The shape of the Eckart potential is shown in Fig. [7J 
We define 

2 2m(E-V± 00 ) T ^oo + ^-oo , Qn ^ 

k ±oo = ^ and k = . (30) 
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The transmission amplitude for this potential is given by (see [27J EE]) 



Y[ika + (1/2) + a/(1/4) - 2mV a 2 /h 2 ] 



x 



r[ifca + (1/2) - v/(l/4) - 2m1/ a 2 /^ 2 ] 



Therefore, the transmission and reflection probabilities are 

sinh(7rfc_ 00 a) sinh(7rfc 00 a) 



(31) 



T 



sinh 2 (vrA;a) + cos 2 (^^(1/4) - 2mV a 2 /h 2 



(32) 



and 



R 



cosh 


na (k - V k l + fc-oo " k 2 ) 


— cos 7r6 


cosh 


vra (k + a/A; 2 ^ + A: 2 ^ - fc 2 j 


— cos 7r6 



(33) 



where k = ^j2/mE~JW and 6 = a/1 — 8mVoa 2 /^ 2 . Fig. [8] shows the transmission and 
reflection probabilities vary as Vo- Fig. [8] (a) and (b) compare the effects of a on the 
probabilities for a = 1 and a = 2 respectively. When a increases, the transmission probability 
approaches unity. Analytically, the transmission resonances occur at 



Va 



h 2 



2ma 2 



n(n + 1). 



(34) 



In contrast, there are no reflection resonances for this potential. 



3.4 A Hulthen potential 

A Hulthen potential takes the form [29] 

V(x) = 6(-x) 



In 



e -ax _ g 



0(x) 



Va 



e ax — q 



(35) 



where Vo, a, and q all are real and positive with q < 1. 9(x) is the Heaviside step function 
defined by 



6{x) 



1, if x > 
0, if x < ' 



(36) 



The shape of the Hulthen potential is shown in Fig. |9j We are interested in two cases 
E > V(x) > and V(x) > E > 0. 
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3.4.1 Case I: E > V(x) > 

The transmission and reflection amplitudes for this potential are given, in terms of a hyper- 
geometric function, by (see [29] ) 

t = 1 q) q [{q(l + 2y){X 2 - 2A/i + /i 2 - is 2 )F(l + A - /i - is, 1 + A - /i + v, 2 - 2/i; q) 

F(X + n - v, A + fj, + is, 1 + 2/i; q)} - {q(l - 2/i)(A 2 + 2A/i + it 2 - is 2 ) 
F(l + A + \x - v, 1 + A + fj, + is, 2 + 2/i; q)F(X - \x - is, A - \x + is, 1 - 2/i; q)} 
-{(1 - 2/i)(2/i)(l + 2/i)F(A + /i - i/, A + /i + i/, 1 + 2/i; q) 
F(X — /x — v, X — /x + is, 1 — 2//; g)}] 

/[{g(A 2 + 2A/i + /i 2 - ^ 2 )F(1 - A - \x - is, 1 - A - \x + is, 2 - 2/i; g) 

F(A - /i - i/, A - /i + i/, 1 - 2//; q)} + {q(X 2 - 2A/i + /i 2 - z/ 2 ) 

F(l + A- /1-Z/, l + A- /i + z/,2-2/i; q)F(-X - \x - v, -A - /i + is, 1 - 2/i; q)} 

-{(2/i)(l - 2/i)F(A - \x - v, A - fx + is, 1 - 2/i; q) 

F(-A-/i-i/,-A-/i + i/,l-2//;g)}] (37) 

and 



- [{(1 + 2/i)F(A + /i - i/, A + /i + i/, 1 + 2/i; g) 



1 + 2/i \ E-k 1 

F(l - X - \x - is, 1 - A - /x + i/, 2 - 2/i; g)} + {(1 - 2/i) 

F(l + A + /i - i/, 1 + A + /i + v, 2 + 2/i; q)F(-X - /x - v, -A - \x + v, 1 - 2/i; g)}] 

/[{g(A 2 + 2A/i + /i 2 - z/ 2 )F(l - A - /i - i/, 1 - A - \x + is, 2 - 2/i; g) 

F{X - \x - v, X - \x + is, 1 - 2/i; g)} + {g(A 2 - 2A/i + ll 2 - is 2 ) 

F(l + X - fx - u,l + X - \x + is, 2 - 2//; q)F(-X - \x - v, -A - \x + v, 1 - 2/i; g)} 

-{(2/i)(l - 2/i)F(A - /x - is, X - fx + is,l - 2/x;q) 

F(-X- Li-is,-X-Li + is,l-2Li;q)}), (38) 

where /i = zfc/a, z/ = zp/a, A = zVo/ag, p 2 = (E + V /q) 2 — m 2 , and k 2 = E 2 — m 2 . The 
transmission and reflection probabilities are derived from 

T = \t\ 2 and i? = |r| 2 . (39) 

We can check that 

T + R = l. (40) 

The transmission and reflection probabilities for varying E are shown in Fig. [TQl There 
are both transmission and reflection resonances. Fig. [10] (a) and [10] (b) describe how the 
diffuseness a has an effect on the probabilities with the other parameters m, Vq, and q fixed. 
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3.4.2 Case II: V(x) > E > 
By the WKB method, we obtain 




Fig. [IT] shows the transmission probability for different heights of the Hulthen potential. 
The probability decreases if the height of the Hulthen potential increases. 

4 Conclusions 

In this paper, we introduced a 2 x 2 transfer matrix for a scattering process (E > V(x)) 
and the WKB method for a tunneling phenomenon (E < V(x)). It is found that the WKB 
method is high accurate when compared to exact solutions. We could use this WKB method 
to find reflection and transmission probabilities for difficult potentials which may give no 
exact solutions. 
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Figure 2: The rectangular barrier with a = 1 and V — 1. 
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Figure 3: The effects of barrier heights Vo of the potential on the probabilities for varying 
q in the rectangular barrier with (a) k = 1, (b) k Q = 2, (c) k Q = 10, (d) k Q = 100 and (e) 
k n = 1000 for o = l. 
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Figure 4: The effects of the widths a of the potential on the probabilities varying as q in the 
rectangular barrier with (a) a = 1, (b) a = 2, (c) a = 10, and (d) a = 100 for k = 1. 
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Figure 5: Plot of the transmission probability and its lower bound. T 2 x 2 denotes the 
lower bound derived from a 2 x 2 transfer matrix. 



15 




Figure 6: The exact solution of the transmission probability compared to the WKB approx- 
imation with (a) V = 1, (b) V = 10, (c) V = 50, and (d) V = 100 for m = H = a = 1. 
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Figure 7: The Eckart potential with — 2, = 1, a — 3, and V = -1/9. 
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Figure 8: The effects of a on the probabilities varying as Vq in the Eckart potential with (a) 
a = 1 and (b) a = 2 for /c_oo = 1, fcoo = 2, and m = h = 1. 
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(a) (b) 

Figure 10: Plotting of transmission and reflection probabilities for varying E for a Hulthen 
potential with (a) a = 0.5 and (b) a = 1 when m = 1, Vq = 1, and q = 0.9. 
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Figure 11: The transmission probability obtained from the WKB approximation with (a) 
Vq = 1, (b) V = 2, (c) Vq = 10, and (d) V = 50 for q = 0.9, a = 0.5, and m = ft = 1. 
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